Brownian motion approach to the ideal gas of relativistic particles 
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The relativistic generalization of a free Brownian motion theory is presented. The global charac- 
teristics of the relaxation are explicitly found for the velocity and momentum (stochastic) kinetics. It 
is shown that the thermal corrections, to the both relaxation times T (of stationary autocorrelations) 
and transient relaxation time of momentum, appear slowing down the processes. The transient re- 
laxation time of the velocity does not depend explicitly on temperature, T(vo) = m(vo)/j = eo/jc 2 , 
and it is proportional to the initial energy of a relativistic Brownian particle. 
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During the last twenty years the Brownian motion ap- 
proach has become the most important tool for discover- 
ing and studying the spectacular phenomena, the critical 
slowing down [1], noise-induced transitions [2], stochas- 
tic resonance [3], ratchet dynamics [4], and resonant ac- 
tivation [5], in (nonlinear) nonequilibrium systems. In 
the most of thermodynamical applications the nonlin- 
ear Langevin equation appears as an overdampcd limit 
of a dissipative dynamic in a certain external potential 
[6]. The important conclusion of the present paper is 
that, within the Markovian diffusion theory, the nonlin- 
ear stochastic kinetic results already on the basic level 
of (spatially homogeneous) ideal gases description if the 
relativistic particles are considered. Particularly the rel- 
ativistic free Brownian motion theory based on identical 
general assumption as the original one [i.e., including the 
interaction with the thermal bath by a systematic force, 
the dissipative kinematic friction proportional to the ve- 
locity — jv, and an additive random force, the thermal 
fluctuation, given by Gaussian white noise £t] 
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where [c — 1 is used throughout the paper] p = 
mv/y/l — v 2 is a relativistic momentum of a Brownian 



particle, 



p/Vp 2 



,2 - 



p/e, e = yV 2 + ' 



and 



(£t£s) = 2D5(t — s), is consistent both with the thermo- 
dynamical requirements of the equilibrium Gibbs ensem- 
ble theory and with the special relativity requirements of 
transformational properties of the momentum probabil- 
ity density distribution. 

First note that the general Kramers-Fokker-Planck 
kinetic equation (see, e.g., the equation (26.26) of the 
Ref. [7]) of a Markovian diffusion description simplifies 
to the form 



d t P = d p L 2 2[(d p Hi)P + k B Td p P] 



(2) 



for the spatially homogeneous system of free noninter- 
acting particles [of the Hamiltonian Hi(p)]. The kinetic 
coefficient L22 = L22 (p, T) is in principle determined 
by the properties of the interaction with the thermal 
bath and do not depend on p for the potential inter- 
action. For the latter case (of the state independent dif- 
fusion) the Langevin description Eq. (1), and the cor- 
responding Fokker-Planck one [8], Eq. (2), are equiv- 
alent in view of the identification L22 = 7 = l{T), 



D = jksT (the Einstein relation, or the fluctuation- 
dissipation theorem), and d p H\ = v. The Eq. (2) is 
the most general form of the diffusion in the momentum 
space leading to the Gibbs-Boltzmann stationary state 
W(p) (x e- 13 " 1 ^ of thermal equilibrium (fi = l/k B T). 
The corresponding (nonlinearly coupled) Langevin equa- 
tion (in the Stratonovich interpretation [8]) is 



Pi 

where L22 = 7ff 2 and g' = dg/dp. 

The (any) corresponding probability density distribu- 
tions W(p) and W'(p') in a resting and moving (with 
a constant velocity V) reference frame are related by 
non-covariant formula eW(p) = e'W'(p'), where e and 
p should be expressed as functions of e' and p' , accord- 
ing to the Lorentz transformation [9]. Because the time 
increments (counted from the initial preparation of the 
system) are given by t' and t 
the relation 
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V 2 , respectively, 



(4) 



Vi -v 2 

is required. One has v = (v' + V)/(l + W), & = (1 — 
V 2 ) 1 / 4 ^', and, using the Stratonovich calculus, dp/dt = 
(1 + Vv')dp' jdt' , thus the transformed stochastic kinetic 
(1) takes the form 

4/ _ ^ v' + V (1 - V 2 fl A 
~d~t i ~ ~ 7 (1 + Vv') 2 + 1 + Vv' 



(5) 



where v' = p' ' / 'y 1 p' 2 + m 2 . It is easy to verify that the 
r.h.s. of (4) solves the Fokker-Planck equation corre- 
sponding to the Eq. (5) if (and only if) P(p, t) satisfies 
the proper Eq. (2) related to the kinetic (1). The same 
conclusion applies also to the general case (3). 

The Eq. (1), p t = — jp/ \/p 2 + m 2 +£ t , for the momen- 
tum or, equivalently, the (Stratonovich) equation for the 
velocity, 

<] 2 fl 2 /m + Z t [l-v 2 fl 2 /m, 
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thus provides the relativistic generalization of the diffu- 
sion Markovian description of equilibration for noninter- 
acting free particles system [10]. The normalized station- 
ary distributions are 

W(v) = - e -(i/D)Vp 2 +™ 2 (j) 

[P > 2mK 1 {m 1 /D) [> 
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and 



W(v) 



(1 - V 2 ) 3/2 ( 7 / D ) m /^l^2 

2K 1 (mj/D) 



(8) 



where 7 = D[3. K n is a modified Bessel function. The 
second moments read 



(P 2 



m Ki(m{3) 



(v 2 ) = 1 



00 **oM 



(9) 



so the (even) moments of the velocity, as incomplete inte- 
grals of Bessel functions, are not given in closed analytical 
form. 

The nonlinear kinetic (1) or (6) docs not belong to 
the class of the known solvable models and the time- 
dependent solution is not known. Nevertheless, there 
are still the quantities, useful to characterize the relax- 
ation of the nonlinear systems, which can be expressed 
by quadratures of the (necessarily autonomous) drift and 
diffusion coefficients, without solving the kinetic equa- 
tion. These are the complete (over time) integrals of 
(various) transient (single-time) moments [11] and sta- 
tionary (two-times) correlation functions [12]. In fact, 
making use of P(x, t\x n ) = e tL 8(x — x ), where L(xq) is 
a backward Fokker-Planck operator [8], and computing 
X(xq) = Jq°° dt(x t ) Xo , where the integrand is a condi- 
tional moment equal xq at the initial time t — 0, by a 
formal integration dte tL = —L^ 1 one obtains 



LX 



-x , 



(10) 



which is a nonhomogeneous linear equation of the first 
order for the function dX/dx - The (x) st = is required 
for the convergence of X and the condition X(0) = 
determines a first constant of integration for Eq. (10). 
The second one should be chosen in such way that the 
correct (deterministic) result is recovered when the noise 
strength D — > 0. The Eq. (10) differs from the Pon- 
tryagin equation for the mean first passage time [13] in 
nonhomogeneous term (— x instead of —1). The inte- 
grals of stationary autocorrelations are then given by 
S x = Jq 00 dt(x t xo) s t — J dxW(x)X(x)x. The last quan- 
tity (i.e., the power spectrum function taken at zero fre- 
quency S x — S x (0)), is called, after some normalization 
(usually divided by the stationary variance) resulting in 
the proper (time) dimension, the relaxation time T [12]. 
The analytical expression for this quantity has been first 
obtained by Jung and Risken [14]. Then, by a simi- 
lar consideration, the formulae for the quantities of the 
former type, called generally nonlinear relaxation times, 
have been found [15, 11]. 

Identifying the appropriate coefficients from Eqs. (1) 
and (6), and carrying out the integrals of the Jung and 
Risken formula [14] we obtain the exact results 
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7 3 Ki(mj/D) 



(11) 
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The same Eq. (11) is obtained for the ordinary Brow- 
nian motion, S° = m 2 S° = Dm 2 /j 2 . Introducing 
T x = S x /(x 2 ) and using m and m/7 as units of mo- 
mentum and time [which is equivalent with considering 
the generic forms of Eqs. (1), (6) with a single (implicit) 
parameter D — r and m = 7 = 1] one gets for dimen- 
sionless quantities 
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where t is a dimensionless temperature, r = D/mj = 
hsT/mc 2 . Note that (13) is the alternative form of 
T v = t/(v 2 ), where (v 2 ) is given by Eq. (9), which is 
more suitable for numerical computations. The respec- 
tive quantities of the ordinary Brownian motion are sim- 
ply equal to the inverse of the relaxation rate of a lin- 
ear model (1), T° = T° = 1, and do not depend ex- 
plicitly on temperature [16]. In contrast both T v or T p 
increases with temperature, see Fig. 1, approaching the 
asymptotic T v (t) s=a r + n/2 or T p (t) (5/2)t, respec- 
tively. The former asymptotic law follows from represen- 
tation (9) [the complete integral of K$ is equal to 7r/2], 
the latter is obtained by the use of the asymptotic for- 
mula K n (l/r) w (n - 1)!(2t)™/2 (for r -> 00) [17] in 
Eq. (14). For small r -> T v (t) w 1 + (3/2)r and 
T p (t) m 1 + 3t. So, the decay of both (stationary) auto- 
correlations in a rclativistic case proceed slower in higher 
temperatures, and slower than for the ordinary Brownian 
motion. We want to stress however that, in the case of 
nonlinear system, the relaxation time T x should not be 
interpreted as a characteristic of the equilibration process 
from a certain nonstationary state. One has in particu- 
lar T x (t = 0) = 1, which would be misunderstood that 
the (deterministic) evolution of the relativistic and clas- 
sical system proceeds with the same (or at least similar) 
speed. Meanwhile it may be generally proved that at the 
limit r — > only the lowest order terms of the drift and 
diffusion coefficients contribute to T x . Thus, such a con- 
clusion applies at most to the relaxation from the close 
to the equilibrium states. 

In order to compare the nonstationary behavior let us 
define the transient relaxation time 



T T (x ) 



r c 
Jo 



dt(x t )a 
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which is related to the quantity X, Eq. (10), via 
x r r (x ) = X T (x ). Then P T (po) satisfies 



PoPt/Jl+rt + TP? = -po. 



The partic ular so lution of the nonhomogeneous equation 
is P' T = y/l +]3q + r, whereas the (general) solution of 
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the homogeneous equation is singular at r = 0. Thus, 

Tr(po) = bo \fl + Vl + log(po + V 1 + Po)]/2po + r. (16) 
Similarly, 

[-«o(l-«o 2 ) 3/2 -3r V o(l-«o 2 ) 2 ]V;' + r(l-«o 2 ) 3 K" = -«b 

has a (regular at r = 0) particular solution V T ' = (1 — 

« 2 )- 3 / 2 , so 

T T (v ) = 1/^1 - vl = s . (17) 

The deterministic (generic) Stokes equations, p = 
—p/ \Jp 2 + 1 and i) = — f(l — v 2 ) 3 / 2 , are solvable by 
quadratures. The solution, with the initial condition xo, 
is 

t{x)=t x (x)-t x {x ), (18) 

where t p (p) = - yj 1 + p 2 + log |(1 + y/l + p 2 )/p|, and 
t„(v) = -(1 - v 2 )- 1 ' 2 + log |(1 + Vl - -y 2 )/w|, respec- 
tively. One can verify that the direct calculation of 
J °° dtx(t)/xo = f" g dxt(x)/xo yields ToO^o) as given by 
Eq. (16) or (17). The previous results (11) and (12) arc 
also recovered when T t (xq) is integrated with XqW(xq), 
where IT is a stationary probability density distribu- 
tion (8) or (7), respectively. Both (dimensionless) tran- 
sient relaxation times for nonrelativistic Brownian mo- 
tion arc equal to unity. The rclativistic relaxation of 
momentum depends both on the temperature and the 
initial state. T T {po), Eq. (16), is an increasing func- 
tion of both arguments and has the following asymptotic 
T T (p Q ) « l+r+^/6, forpo -> 0, andT r (p ) « M/2+r, 
for po — * oo. The result for the velocity is somehow un- 
expected. The Eq. (17) shows that indeed the relativistic 
relaxation proceeds slower, however does not depend on 
temperature. It turns out that with the increase of tem- 
perature the early stage of evolution becomes faster in 
such a way that compensates the subsequent slower long- 
time behavior, see Fig. 3. The transient relaxation time 
T t (vq) is equal to the (dimensionless) initial energy of 
the relativistic Brownian particle The plots of To(xo) 
are shown in Fig. 2. 

The equation (1) (in a generic form) has been solved 
numerically by the Runge-Kutta method for 10 4 sample 
realizations of the white Gaussian noise, for different ini- 
tial conditions and temperatures. (p t ) and (v t ) have then 
been computed as the arithmetical average of appropri- 
ate values obtained for different trajectories of the noise. 
The typical curves are plotted in Figs. 3 and 4, together 
with the corresponding deterministic results (18). The 
solution of the relativistic Stokes equation for velocity 
changes qualitatively, when vo exceeds 1/2. In fact be- 
low this value of the initial state the function v(t) is con- 
vex down for all t [as well as (always) p(t)]. However for 
vq > 1/2 at the early stage of evolution the v{t) appears 
convex up [in order to achieve the required asymptote 



v(t) = 1 = c for vo — > 1]. The former subrelativistic case, 
with particular v = 2/5 (po ~ 0.436), and the latter, 
with po = 4 (v ~ 0.970), are presented. 

The results of the paper may be summarized as fol- 
lows. We have shown that using the relations of the rel- 
ativistic dynamics we can generalize the Brownian mo- 
tion theory in a way consistent with thermodynamical 
requirements. The general solvable equation (10) for the 
quantities simply related to the transient relaxation times 
(15) has been found. These, as well the stationary re- 
laxation times T, have been found analytically for the 
rclativistic momentum and velocity (stochastic) kinetics, 
Eqs. (13, 14, 16, 17). All these quantities (considered 
as dimensionless) are equal unity for the ordinary Brow- 
nian motion. The relaxation times T x (t), by the nor- 
malization equal unity in the deterministic limit r — > 0, 
exhibit positive first order corrections for small temper- 
atures, becoming just proportional to t in a high tem- 
perature asymptotic. Thus the decay of correlations in 
a stationary state proceed slower than for the nonrela- 
tivistic Brownian motion. The transient relaxation time 
[describing the relaxation from a given initial state xq, 
and by the normalization equal unity for the determin- 
istic system set initially at the equilibrium, Tb(0) = 1] 
of momentum T T (po) is an affinic function of tempera- 
ture. The asymptotic behavior for the large argument(s) 
reads T T (p ) « \Po\/2 + T - The transient relaxation time 
for the velocity (17) does not depend on temperature, 
and it is equal to the dimensionless initial energy of the 
relativistic Brownian particle. The thermal corrections 
to the considered global characteristics of the relaxation 
are linear with the slope of the unity order and appear 
(even) for the noninteracting free relativistic Brownian 
particles system. The corrections can be observed for suf- 
ficiently high physical temperatures. For, e.g., hydrogen 
particle (proton) the dimensionless quantity, which char- 
acterize the deviation from the ordinary nonrelativistic 
behavior, r = k B T/mc 2 w 10" 13 T[K]. So for the tem- 
perature T ~ lMeV ~ 10 10 K the relative correction is 
of the order 0.1%. The correction to transient relaxation 
time of the velocity has the kinematic origin only and 
for T T (po) both kinds of corrections are decoupled. Note 
that the general dependence on temperature in physi- 
cal units, introduced (multiplicatively) by the coefficient 
7(T), should be counted separately [16]. 

The three-dimensional generalization of the proposed 
diffusion description is obvious, however in contrast to 
the ordinary Brownian motion theory it leads to the set 
of the three coupled (v = p/e = pj y/ m 2 + p 2 ) Langevin 
equations (with statistically independent noises in each 
direction) or equivalently to the (forward) Fokker-Planck 
description [8] [V = d/dp\ 

8 t P(p, t\p- Q ) - 7 V o (pP/e) + DV 2 P. (19) 

The required form of the (Gibbs-Boltzmann) stationary 
distribution is still recovered, however the calculations of 
relaxation times cannot be analytically completed in this 
case. 
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FIG. 1: Plots of T p [Eq. (14), upper curve] and T v [Eq. (13)] 
vs t (dimensionless temperature). The latter has been com- 
puted numerically by the Romberg method. 
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FIG. 2: Plots of T (po) [Eq. (16), r = 0, upper curve] vs p 
and T T (v ) [Eq. (17)] vs 10xw . 



FIG. 3: Plots of {vt) VQ /vo vs t. The upper and lower curves 
correspond to vo ~ 0.970 and vo = 2/5; respectively. The de- 
terministic solutions, Eq. (18), are marked by (r = 0). The 
dotted curve, e~', corresponds to the nonrelativistic universal 
result. The remaining curves have been obtained numerically 
for t = 0.2 and 0.6, respectively. 



FIG. 4: The same for {pt) PO /po vs t. 
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